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.ABSTRACT 


An  approximation  to  the  Debye  characteristic  temperature  at  0  °K 
is  developed  and  applied  to  several  cubic  metals.  The  results  obtained 
are  estimated  to  be  accurate  to  better  than  1$# 

In  addition,  the  Bom  and  Begbie  theory  of  lattice  vibrations  is 
applied  to  the  calculation  of  the  specific  heat  and  Debye  temperature 
of  silver  over  the  temperature  range  0°  -  150°K.  Three  such  calculations 
were  carried  out;  the  Houston  approximation  to  the  frequency  spectrum 
used  in  the  first  calculation  is  extended  in  the  second  calculation  to 
take  into  account  the  variation  with  temperature  of  the  elastic 
constants, and  is  modified  in  the  third  calculation  so  as  to  give  the 
correct  results  at  both  low  and  high  temperatures.  The  results  agree 
with  experiment  to  approximately  the  same  degree  as  the  similar 
calculations  of  Bhatia  and  Horton  on  the  basis  of  the  Bhatia  theory. 
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INTRODUCTION 


The  subject  matter  divides  naturally  into  three  fairly  well- 

defined  regions  of  discussion  and  is  so  treated.  In  Part  I*  after  an 

outline  of  the  historical  development  of  specific  heat  theory,  a 

general  expression  for  specific  heats  is  derived  and  is  discussed  at 

some  length,  An  accurate  approximation  to  the  Debye  characteristic 

temperature,  ©5  ,  of  cubic  metals  at  0°K  is  developed  in  Part  II  and 

is  applied  to  the  calculation  of  0D  of  several  common  metals.  In 

Part  III  Born  and  Begbie ' s  theory  ©f  lattice  vibrations^  is  applied 

to  the  calculation  of  the  specific  heat  of  silver.  Three  calculations 

of  this  type  are  carried  out,  with  an  attempt  in  two  of  them  to  take 

into  consideration  the  variation  with  temperature  of  the  elastic 

constants.  The  results  are  compared  with  experiment  and  with  similar 

2 

calculations  by  Bhatia  and  Horton  based  on  the  Bhatia  theory  of 
lattice  -coftducl-fdnelectron  interactions  in  monovalent  cubic  metals^. 
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PART  I.  HISTORICAL  DEVELOPMENT  OF  SPECIFIC  HEAT  THEORY 
Classical  Theory 

On  the  basis  of  classical  statistical  mechanics,  the  energy  of 
an  oscillator  moving  under  harmonic,  i.e«  linear,  forces  is  KT  per 
degree  of  freedom,  where  K  is  the  Boltzmann  constant  and  T  is  the 
absolute  temperature*  Since  there  are  three  degrees  of  freedom  for 
each  particle,  the  internal  energy  U  of  a  system  of  N  particles 
is  given  by 

(1)  U  -  3  N  K  T  . 

The  heat  capacity  at  constant  volume ,  C^  ,  is  obtained  by  taking  the 
derivative  of  the  internal  energy  with  respect  to  temperature*  Thus 

(2)  CT  -  3  N  K  . 

This  is  in  reasonably  good  agreement  with  experiment  at  high 
temperatures  but  fails  completely  at  low  temperatures  where  Cv  falls 
to  low  values*  The  introduction  of  quantum  theory  corrects  this 
situation* 

Ehergy  and  Heat  Capacity  of  a  System  of  Harmonic  Oscillators  -  Quantum 

Theory 

An  elementary  result  of  quantum  theory  is  that  the  energy  of  a 
one-dimensional  harmonic  oscillator  is  given  by 

(3)  Wn  -  (n  +  i)  hw 

where  h  is  Planck  *s  constant  divided  by  2tt, 

co  is  the  angular  frequency  of  vibration,  and 
h  is  the  quantum  state  of  the  oscillator* 

The  probability  that  a  given  oscillator  in  thermal  equilibrium  will 
be  in  state  n  (energy  state  W  )  is 
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average  energy  Wn  of  such  an  oscillator  is 
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For  high  temperatures  (hco  «  KT),  the  expression  reduces  to  the 
classical  result  Wn  *  KT. 

The  internal  energy  of  a  system  of  N  such  oscillators ,  each 
with  three  degrees  of  freedom  is 

(6)  U  -  3N  Wn 

*  (3/2)hto  +  3N  hu>  / (e^00^'1  -1) 

Differentiation  of  (6)  with  respect  to  T  gives  the  heat  capacity 

(7)  c  -  3H  (R^)2-^"/KT 

KT2  (eh“/KT  -l)2 

-  3NK  x2ex  /(ex  -l)2 

-  3NK  E(x) 

where  x  =  hco/KT  ,  and 

E(x)  =  x2ex  /(ex-l)2  is  the  Einstein  function. 
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It  was  proposed  "by  Einstein  to  use  this  expression  for  as  an 
approximation  to  the  entire  heat  capacity  of  a  solid.  At  high  tempera¬ 
tures  (x— >0),  E(x)  approaches  1  and  goes  to  3HK,  the  classical  result. 

At  low  temperatures  E(x) ,  and  hence  Cv,  fall  to  low  values .  Though  the  result 
is  qualitatively  correct,  it  fails  to  give  correctly  the  details  of  the 
approach  of  C^.to  zero. 

Frequency  Distribution  Function 

It  is  necessary  to  replace  the  single  frequency  used  in  the  Einstein 
theory  by  a  range  of  frequencies,  characteristic  of  the  substance,  whose 
individual  contributions  to  the  specific  heat  are  expressed  in  a  freqency 
distribution  function  G-(u))  do),  the  number  of  frequencies  lying  between 
a)  and  a)  +  dot).  The  internal  energy  is  then  given  by 

(8)  U  =  |  &(»)to(eWKT  _  l)-idu) 

with  the  condition  that 

(9)  7  |  &(»)do)  =  3U0V 

where  is  the  total  number  of  modes  of  vibration  in  volume  V,  and 

is  the  number  of  atoms  per  unit  volume. 

Differentiation  of  (8)  with  respect  to  T  leads  to  the  general 
expression  for  the  specific  heat  c^  per  atom  at  constant  volume 

(10)  H  7c  =  Kfe/KT)2  f  G(u>)<D2eWKT  (ehio/KT  _  i)-2d;u 

o  v  J 

Determination  of  G-Qp) 

Consider  a  cube  of  side  L.  The  condition  that  a  vibration  wave 
be  periodic  within  the  cube  is  that  there  be  an  integral  number 
of  wavelengths  X  between  the  ends  of  the  cube.  Since  k  =  2rr/X,  the 
necessary  condition  for  periodicity  becomes  that  k^L,  k^L,  and  kzL  be 
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multiples  of  2rr.  The  allowed  values  of  the  wave-' vector  k  may  he 
represented  hy  the  points  of  a  simple  cubic  lattice  in  k- space  with 
lattice  constant  2rr/L.  The  number  of  points  in  an  element  of  k-spaee 
dkxdkydkz  is 


(id 


dkxdky  dk  ^ 
(2rr/L)3 


(2tt)3 


dkxdkydkz . 


In  polar  co-ordinates,  the  element  of  volume  in  k- space  is  given  by 
k2dkda,  where  do.  is  an  element  of  solid  angle*  The  number  of  modes 
between  k  and  k  +  dk  in  a  solid  angle  da.,  N(k)dkda,  is 

(12)  H  (k)  dkda  =  Y  k2dkdfL 

(2n)^ 

=  — ijc 2®£.  dooda 
(2ny  da) 

=  g®  (a))  diada. 

since  k  is  a  function  of  a>.  Waves  of  a  given  k  occur  with  three 
different  polarisations  (one  longitudinal  and  two  transverse),  so  that 


3 


The  required  frequency  distribution  function  G-OiOdcodo.  is  obtained  by 
integrating  (13)  over  all  solid  angles,  giving 

(l^)  VGr(a>)dco  =  ■  dco  f  y  k^(a). )  )  dil  . 

(2rr)  J  tr  1  to>± 

Thus,  on  the  basis  of  a  given  model,  G-(a))  can  be  determined  in  principle 
from  the  dispersion  (00  vs,  k)  relations. 
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Debye  Approximation 

If  the  solid  is  regarded  as  continuous,  then 


(15) 


>i  =  nk 


where  the  velocity  of  the  wave  Vj_  is,  in  general,  a  function  of  direction. 
G-(uj)  becomes 

(l6)  G(u))  da)  =  (l/2rr)^o>2da )  \ Z  da/v3(9,(P) 

i-\  1 


with 

(17) 

Since 


c<1  =  (l/2n)3  Jz  d 
G(u>)du>  =  cq  J 


u)  du>  would  he  infinite  if  to  were  allowed  to  assume 


all  values,  we  require  that  there  he  a  cut-off  frequency  «)m  such  that 


(18) 


(x(co)do)  =  OCjjD^djw  for 


U)  =  O), 


=  0 


for  a)  >. 


'm  * 


m* 


I  Gr(o))du)  =  3N0. 


The  maximum  frequency  is  fixed  hy  the  requirement  that 
Therefore 

(19)  oum  =  (9ti 0/<X1)1/3. 

Introducing  (l6)  for  S(oi)  in  (10),  replacing  0<^  by  9H0/®^  and  i»/K5  by 
x,  and  rearranging,  one  obtains 


(20)  N0Vcr  = 


r 

TO(Kl/fm)3  J0  x2U(x)dx. 

Since  hu)/K  has  the  dimensions  of  temperature,  we  can  write 

(21)  eD  =  tam/K. 

This  is  the  definition  of  the  Debye  characteristic  temperature.  (20), 
written  in  terms  of  9^,  becomes 

(22)  cv  =  9K(T/eD)3D(9I)/T) 
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where 

rO p/T  p 

(23)  D(9D/T)  =  J0  x^Etodx. 


Low  Temperature  Behaviour  of  Debye  Approximation  to  cv 

As  T  ->  0,  Op/T  l,and  the  upper  limit  of  integration  in  D(0^/T) 
may  be  taken  as  infinity. 


(24)  D(eD/T)T^0 
Consider  the  integral 

(25) 


=  If  xV(ex  -  l)-2d*. 


=  r^-D-i 


=  1 


4  e  -  1 


co  r 

>  1 

o  4  J 


dx 

00 


xV 


dx 


(eA  -  1): 


=  L  D  (6jT). 
4  u 


Now 

(26) 


dx 


e-^dx 


=  t»4/15  . 

Therefore  D(<x>)  =  WVl5 
and 


(27)  c  =  9K(T/©J3  4tt^/15. 

v  5^0  D 

Thus  at  low  temperatures,  cv  obeys  a  T^  law. 


Failure  of  Debye  Approximation 

The  range  of  validity  of  (27)  is,  unfortunately,  not  too  great. 

The  limit  on  the  usefulness  is  set,  not  by  mathematics,  but  by  physical 
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reality.  The  Bebye  approximation  just  does  not  represent  the  true 
picture.  For  instance,  since  0^  =  tlu)m/Z  =  (^/^J^^h/E,  ©d  should 
he  independent  of  temperature*  However,  if  one  substitutes  actual  c  1 s 
into  (22),  one  obtains  a  0^  which  is  not  constant  with  temperature,  but 
exhibits  a  minimum.  The  reason  for  the  failure  of  the  Debye  theory  is, 
of  course,  due  to  the  neglect  of  the  discrete  structure  of  matter.  Since 
the  advent  of  the  Debye  theory,  more  elaborate  treatments  of  specific 
heats  have  arisen,  which  take  into  account  the  atomic  nature  of  matter. 


Born  -  von  Karman  Theory 

Until  recently,  it  has  been  customary  to  determine  the  frequency 

distribution  on  the  basis  of  the  Born  -  von  Karman  model*5,  assuming 

central  Hooke*  s  Law  forces  between  an  atom  and  its  nearest  and  next 

nearest  neighbours.  This  introduces  two  unknown  force  constants  which 

are  then  determined  by  the  requirement  that  the  velocities  of  long 

elastic  waves  as  determined  from  elasticity  theory,  and  expressed  in 

terms  of  the  elastic  constants,  correspond  to  the  velocities  as  obtained 

from  the  microscopic  theory.  There  are,  in  general,  three  independent 

elastic  constants  for  cubic  crystals  (cpp,  cpg*  ani  C44) .  One  relation 

between  the  elastic  constants  exists  if  either  the  crystal  is  elastically 

isotropic  (2c ^  =  c^  -  c1  ,  or  the  entire  interactions  a,re  central  and 

the  Cauchy  relation  is  satisfied  (c.^  =  .  Unfortunately,  in  metals 

this  is  seldom  the  case.  Hence,  one  may  use  an  approximation, 

6 

introduced  by  Fine  ,  by  means  of  which  the  three  elastic  constants  are  reduced 


to  two  independent  ones. 
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Bhatla  Theory 

More  recently,  Bhatia^  has  developed  a  new  theory,  taking  the  cotduchon 


electrons  into  account  in  an  approximate  way.  This  introduces  a  non¬ 
central  component  into  the  interactions  and  yields  three  atomic  force 
constants  which  can  he  uniquely  related  to  the  three  independent 
elastic  constants. 

Born  and  Begbie  Theory 

In  19 47,  Bom  and  Begbie^  extended  the  Born  -  von  Karman  model  to 
include  non-central  forces  (general  forces)  and  it  is  this  theory  which 
will  he  applied  here  to  the  calculation  of  the  specific  heat  of  silver. 
If  one  applies  the  Born  and  Beghie  theory  to  face-centered  cubic 
crystals,  such  as  silver,  one  finds  that  the  nearest  neighbours 
introduce  three  force  constants  which  can  he  uniquely  determined  in 
terms  of  the  three  elastic  constants. 

Secular  Bquation  for  the  Frequencies 

On  the  basis  of  the  Born  and  Beghie  theory  of  lattice  vibrations, 
the  secular  equation  for  the  frequencies  is 


Di  j  = 

I±5  =1  for  i  =  3, 

=  0  for  i  f  j, 

and  =  cosCak^/ /2)cos(aky/2) 

=  sin(aki//2). 


*  a1  is  the  nearest  neighbour  distance  and  the  force  constants  and  elastic 
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constants  are  related  in  the  following  way: 

(29)  <*•  =  0,44  “  °li/2 

ft  =  cxi^2 

V  =  ^cxx  —  °44  ^  )/2 

with  6  =  ou  -  c12  -  Z%k 

o 

The  secular  equation  (28)  has  three  solutions  of  for  a  given 
k  direction,  corresponding  to  one  longitudinal  end  two  transverse 
“branches.  However,  we  should  note  that  only  in  certain  special 
directions  are  the  solutions  purely  longitudinal  and  transverse. 


.  ;  '  £ 
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PART  II.  DEBYE  Qd  AT  LOW  TEvQPERATURIS 

Validity  of  Debye  Approximation  at  Low  Temperatures 

We  have  found,  using  the  Debye  approximation,  that 

(i)  eD  =  (sNoMiJ^Vk. 

This  expression  can  be  used  to  calculate  ©p  temperatures  near  0°K, 
since  only  the  lowest  frequencies,  or  equivalently,  the  longest  wave¬ 
lengths,  are  excited  in  this  temperature  region.  The  lattice  structure 
of  the  metal  is  not  apparent  to  long  waves,  so  that  a  continuous 
medium  approximation  is  applicable. 

Approximations  to  o(^ 

We  now  proceed  with  the  evaluation  of  which  is  given  by  ( 1.17 ;  as 


Since  crystals  in  general  are  not  isotropic,  the  vi  depend  on  direction 
in  a  complicated  way,  with  the  result  that  the  integral  must  be 
evaluated  by  approximate  methods.  Several  such  methods  have  been  devised 
and  we  shall  discuss  seme  of  them  in  greater  or  lesser  detail. 

Hopf  -  Lechner  Approximation 

Until  recently,  a  method  often  used  has  been  that  of  Hopf  and 
7 

Lechner  .  The  integral  (2)  is  evaluated  by  expanding  the  integrand  in 
a  polynomial  in  vf  and  determining  the  coefficients  by  the  Lagrange 
interpolation  method.  Two  of  the  defects  of  this  approach  are  that  the 
method  is  extremely  labour ious  and  that  there  is  no  means  of  ascertaining 
the  accuracy  of  the  result. 
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Q,uimby  -  Sutton  Approximation 

A  modification  of  the  Hopf  -  Lechner  method  is  that  of  Quimby  and 
Sutton8.  It  is  based  upon  an  evaluation  of  the  distribution  in  solid 
angle  of  the  roots  of  the  cubic  equation  which  describes  the  propagation 
of  acoustic  waves  in  the  crystal.  The  authors  claim  an  uncertainty  in 
the  computed  value  of  %)  due  to  the  approximation  of  the  order  of  0.1$ 
when  the  anisotropy  factor  ^  =  2044/(0x1  -  042)  is  as  large  as  4.  The 
method  does  appear  to  involve  a  considerable  amount  of  work. 

Blackman  Approximation 

A  very  simple  semi-theoretical  formula  for  %),  designed  to  cover 
cases  where  ^  is  considerably  greater  than  1,  has  been  developed  by 
Blackman9.  Calculations  can  be  quickly  made  with  the  formula  to  an 
accuracy  of  5  -  10$. 

Bhatia  -  Tauber  Approximation 

More  recently,  a  simple  approximate  expression  for  the  integral  (2) 
was  derived  by  Bhatia  and  Tauber10,  based  on  a  method  employed  by  Houston11 
to  obtain  the  approximate  form  for  the  frequency  spectrum  of  crystals. 

The  integrand  is  expanded  in  those  Kubic  Harmonics  first  introduced 
by  von  der  Lage  and  Bethe12,  which  are  invariant  under  the  operations  of 
the  cubic  symmetry  group.  Thus 


(3) 

<*1  =  u/2*)3 

Since 

(4) 

J 

K4*  =  0 

m  £  0 

*  4tt 


m 


0 


)  i.  >  -L1,  ■  "■'/  ...  ; . 


i'.  ;  o  w.  h : 

vb:>  ;  1  Ui.  $0  -  cQ!  ■ :  r  ii*  ‘3  'j  0  "  ;!./ 

■  o  \  i  . 
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we  have 

(5)  »  (1/2tt)°4it^0. 

If  the  crystal  is  not  too  anisotropic,  the  first  three  Kubic  Harmonics 
are  sufficient  to  give  a  good  approximation  to  the  integral.  In  terms  of 
the  direction  cosines  x,  y,  z, 

(6)  K0  =  1 

*  x4  +  y4  -t-  z4  -  3/5 
Kg  -  sW  +  (1/22)^  -  1/105 

The  coefficients  j30t  and  fio  are  determined  analytically  from  the 
solutions  of  the  secular  equation  (1.28)  along  the  three  principal 
directions  of  a  cubic  crystal, (100) ,  (110),  and  (111).  The  quantities 
referring  to  these  three  directions  will  be  distinguished  by  the 
subscripts  A,  3,  and  0,  respectively.  is  then  given  by 
(V)  /30  *  (1/35)  H  {lO(l/vf)A+  16(1/v®)b  +  9(l/vf)0  }. 

C-f 

Therefore, 

(8)  ^  =  (l/2ir)3(4n/35)  £  (lO(l/v^)A+  16(l/vf)B+  9(l/vf)0]  . 

{-! 

A  better  approximation  to  ^  can  be  obtained  by  using  more  terms 
in  the  expansion  in  Kubic  Harmonics.  By  taking  three  more  directions, 
(210),  (211),  and  (221),  denoted  by  3,  E,  and  E,  respectively,  Betts, 
Bhatia,  and  Wyman^  have  derived  the  following  six-term  approximation 
to  d  1  : 

1  3 

(9)  =■  (l/Sir)3(4ir/l,081,080)  £  [ll7,603(l/v3)A  +  76,544{lAf  )B 

C-f 

*  17,496(l/vJ)c  *  381,350(l/vf)D  +  311,040(l/vf ) B  + 
177,l«(l/^)s.}. 

A  more  complete  discussion  on  Kubic  Harmonics  and  their  properties  will 
be  found  in  Ref.  13. 


_  - 
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Schiff  Approximation 

At  the  suggestion  of  Dr.  Schiff ,  the  above  approximations  may  he 
improved  up  on  in  the  following  way. 

On  expanding  the  determinant  (1.28),  one  obtains  a  cubic  in  cw2 
(10)  a)6  +  aa/^  +  ba>2  +  c  =  0 

where  a,  b,  and  c  are  functions  of  the  S^,  and  the  elastic 

constants.  The  solutions  are  uj?(k ,©,<#>)  where  i  =  1,  2,  3*  In  the 

2 

long  wave  limit  =  v^k,  so  that  the  cubic  equation  in  a)  may  be 
2 

written  as  one  in  v  ,  namely 
(Xl)  V6  -  Lv^  +  Mv2  -  N  =  0 

with  solutions  v2,  v2,  and  v2  which  are  functions  of  direction.  In 
terms  of  the  solutions  of  (ll) ,  the  coefficients  are  given  by 
(12)  L  =  v2  +  v|  +  v^ 

M  =  v^v2  +  v|v^  +  v^v| 

H  =  v2v|v2  . 

Now  consider  the  required  integral  (2) 

(/I  +  I  +I  W=  f  (^3) 3  +  (^j)3  ♦  ^2)3  a*. 

-H  V1  Vf  J  <V2T3)3 

It  can  be  rewritten  as 


(13) 


where 


"  ^V2V3^  +  ^v3vl^3  +  (vlv2^  ((vgT^)2*^.^)2  +  (vjVg)2]^  2 

[(v2v3)2  +  (^1)2+  (^)2]3/2  [(t2v2t2)]3/2 

=  Je(M/U)3/2cL  JZ 


q  q  q 

r  =  (v2v3^J  +  »  (viv2' 

[(V2V3)2  +  (v^)2  +  (v-jVg)  2]3/2 
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It  turns  out  that  R  is  a  very  slowly- varying  function  of  angle, 

(R  must  lie  between  1  and  1//3)  •  Por  example,  for  the  worst  case  in 
a  fee  metal,  namely  Pb,  R  lies  between  0.62  and  0.76.  Because  of  the 
small  speead  in  R  over  the  unit  sphere,  we  expect  it  should  be  possible 
to  obtain  a  good  approximation  to  by  writing  it  as 

(15)  «l1  =  (1/210%  J(M/H)3/2a m. 

where  £  is  some  appropriate  average  value  of  R  over  the  unit  sphere. 

The  exact  expression  for  5  is  of  course  given  by 


(16) 


jR(M/iir)3/2a.a 

J  (M/ir)3/2dj2. 


The  numerator  in  the  expression  is  just  the  integral  we  are  trying  to 
get.  However,  we  now  find  an  approximate  value  for  S  by  expanding  the 
integrands  of  both  numerator  and  denominator  in  Rubio  Harmonics,  using 
only  the  first  three  terms.  Thus 


(17) 


approx. 


10%(M/H)j/2  +  16%(M/H)g/2  +  9$c(H/H)g/2 

10(M/H)3/2  +  16(M/IT)3/2  +  9(M/H)3/2 


=  p/q. 

Since  R  has  a  small  spread  over  the  unit  sphere,  P  and  Q,  vary  in 
approximately  the  same  way,  and  the  errors  in  P  and  Q,  introduced  by 
the  approximation  are  about  equal.  Thus 

(18)  E  =  ill  ±±1  (6,~et) 

Q(1  +  fjJ 

ss|  (1  +  *,)(l  -  *t) 

4*  |  (l  +  «,  - 

How  P  is  nothing  more  than  the  Bhatig,  -  Tauber  (B  -  T)  formula  for  ot  , 
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and  is  the  error  introduced  by  that  approximation;  hence,  provided 
that  the  integral  in  (15)  can  he  accurately  evaluated,  we  would  expect  this 
approximation  to  he  considerably  better  than  the  three-term  B  -  T 
approximation  for  cLi ,  The  accuracy  of  this  method  may  he  extended  even 
further  by  now  using  the  six- term  expansion  for  P  and  Q. 

Evaluation  of  S 

The  task  of  evaluating  U  is  made  quite  simple  if  one  first 
calculates  the  required  velocities,  then  Inserts  them  into  the  expression 
for  E,  rather  than  expressing  the  numerator  and  denominator  of  E  in 
terms  of  the  elastic  constants.  The  complexity  of  the  expression 
derived  on  the  basis  of  the  latter  approach  is  easily  appreciated  if 
one  but  looks  at  the  general  formulas  for  the  velocities  in  terms  of 
the  elastic  constants  in  the  directions  D,  1,  and  E. 

Data  Used  in  the  Calculations 

The  constants  used  in  the  calculations  are  given  below.  The 
sources  of  the  data  are  included  as  well,  and  in  the  case  of  the  elastic 
constants,  the  temperature  at  which  they  were  measured.  The  metals  are 
arranged  in  order  of  increasing  anisotropy. 
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table  1. 


Metal 

Density  p 

Elastic  Constants 

H  —  A vo  gadro  *  s  Ho  ^ 
At.Wt. 

Value 

Source 

°n 

C12 

c44 

Temp. 

Source 

Value 

Source  of 
At.  Wt. 

W 

19.0 

HOP 

50.1 

19.8 

15.1 

RT 

Kittel  14 

6. 222 

HCP 

A1 

2.699 

HCP 

12.30 

7.08 

3.090 

o°k 

Sutton  15 

6.030 

HCP 

Ee 

7.88 

HOP 

23.7 

14. 1 

11.6 

RT 

Kittel  14 

8.502 

HCP 

Ag 

10.492 

HCP 

12.0 

8.97 

4.37 

29l°K 

Roehl  16 

8 

00 

6 

HCP 

Cu 

8.9326 

HCP 

17.0 

12.3 

7.52 

RT 

Kittel  14 

8.^70 

HCP 

Pb 

11.34 

HCP 

4.8 

4.1 

1.4 

RT 

Kittel  14 

3.297 

HCP 

K 

0.87 

HCP 

0.46 

0.37 

0.26 

RT 

Kittel  14 

1.341 

HCP 

Ha 

1.0066 

HCP 

0.970 

0.830 

0.58O 

o°k 

Euchs  17 

2.6372 

HCP 

Remarks:  The  density  is  in  gm/cc,  the  elastic  constants  are  in  units 

1 1  22 
of  101  dynes /cm2,  and  H  ,  the  number  of  atoms/cc,  is  in  units  of  10  / cc. 

HOP  is  the  "Handbook  of  Chemistry  and  Physics,  35tk  Edition  (l953~54) ”• 

RT  stands  for  room  temperature,  and  At.  Wt.  for  atomic  weight. 
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General  Expressions  for  the  Velocities  of  Long  Waves 

For  long  waves  ka //2  is  very  small,  so  that  in  the  secular  equation 

(1.28)  we  can  replace 
sin  ka//U  by  ka//2 

cos  ka//2  by  1  -  k2a2/4  , 

Making  these  substitutions  and  writing  v~  for  /kw  in  the  solutions  of 

(1.28) ,  we  obtain  the  velocities  of  long  elastic  waves  in  a  culbic  crystal. 

(19a)  Direction. A  (100) 

/JT1  =  ®11 
l°Ytl  -  °44 
f>7tS  *  °44 

(19b)  Direction  B  (110) 

/>71  "  cll-^ 

/°7tl  =  °44 
/>'rt2  ~  C44  *  s'  * 

(19c)  Direction  G  (111) 

f>\  “  °ll  “  (2/3)  k 

r\x  -  °44  *  (1/3)  fe 
/°Vt2  =  O44  -  (1/3)  * 

Direction  D  (210) 

/°V1  1  °ll/2  ♦  °44/2  +  {/9(°ll  -  C44)2  t'16(°12  +  °44)a  )/10 
f  vtl  *  °44 

r™  "  ,V*  ^W2  -  <Ml  “  °44)a  ^16(c12^  °44)a  1/10 


(194) 
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(19a) 

Direction  ] 

3  (211) 

rl  - 

C44 

*  (5°lf 

°12_4044)/12  + 

(i(3° lX-ei2 

■4044) 

2 +  32(0.  + 

°44) 

S)/12 

rh- 

°44 

+  6/6 

rlz  - 

°44 

*  (5011+ 

C12-4C44)/12  - 

-4044) 

2+32(012+ 

°44} 

S)/12 

(19f) 

Direction  F  (221) 

rl  - 

C44 

*  (5°11+ 

40ia-°44)/18  + 

(/(3°  ll+4012+044 

)2+32(oi;j  + 

°44> 

2)/b 

rh  - 

°44 

+  (5°lF 

4012-°44)/18  - 

(/(3°ll+4°l2*C44 

>a+32(°12+ 

°44> 

2)/l3 

rlz  ■- 

C44 

+  4«/9 

Velocities  of  juong  Elastic  Waves  in  Several  Cubic  Crystals 

The  velocities  for  W,  Al,  Fe,  &g,  Ou,  Pb,  K,  and  Na,  as '  calculated 
from  ^19)  using  the  constants  of  Table  1  are:  (in  units  of  10  cm/ sec) 

TABLE  2. 


(a)  Tungsten 


Dir. 

vi 

vtl 

H2 

100 

5.1350 

2.8191 

2.8191 

no 

5.1325 

2.8191 

2.8238 

ni 

5.1316 

2.8222 

2.8222 

{e; 

Iron 

Dir. 

vi 

vti 

vt2 

100 

5.4842 

3o8368 

3.8368 

no 

6.2214 

3.8368 

2.4681 

in 

6.4484 

2.9947 

2.  9947 

210 

6.0195 

3.8368 

2.9263 

211 

6.2651 

3.4415 

2.9020 

221 

6.3576 

3.4710 

2.6552 

(b)  Aluminum 


Dir. 

'vi 

vti 

vt2 

100 

6. 7507 

3 . 38o 7 

3.3837 

no 

6.8812 

3.3837 

3.1097 

in 

6.9241 

3.2036 

3.2036 

(d)  Silver 


Dir. 

Vtl 

vt2 

100 

3*3820 

2.0409 

2.0409 

no 

3.7628 

2.0409 

1.2017 

in 

3.8814 

1.5333 

1.5333 

£lo 

3.6517 

2.0409 

1.5057 

211 

3.7810 

1.8050 

1.4877 

221 

3.8329 

1.8249 

1.3215 

v'Xtfe 
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(e)  Copper 


Dir. 

V1  vtl  vt2 

100 

4.3625  2.9015  2.9015 

110 

4.9819  2.9015  1.6220 

111 

5. 1719  2. 1354  2. 1354 

210 

4.8119  2.9015  2.0726 

211 

5.0180  2.5474  2.0492 

221 

5.0958  2.5744  1.8094 

(g)  Potassium 


Dir. 

V1  VU  Vt2 

100 

2.2994  1.7287  1.7287 

110 

2.7854  1.7287  0.71919 

111 

2.9296  1.1580  1.1580 

210 

2.6675  1.7287  1.0772 

211 

2.8220  1.4713  1.0658 

221 

2.8742  1.4872  0.88985 

(f)  Lead 


Dir. 

V1  Vtl  vt2 

100 

2.0574  1.1111  1.1111 

110 

2.2713  1.1111  O.55555 

111 

2.3383  0.78567  0.78567 

210 

2.2073  l.llll  0.77152 

211 

2.2804  0.96225  0.75872 

221 

2.3106  0.97554  0.64152 

(h)  So  dim 


Dir. 

V1  Vtl  *t2 

100 

3.1043  2.1004  2.4004 

110 

3.8344  2.4004  0.83391 

111 

4.0487  1.5441  1.5441 

210 

3.6635  2.4004  1.4060 

211 

3.8923  2.0182  1.3919 

221 

3.9673  2.0402  1.1218 

With  the  v.  calculated,  it  is  now  hut  a  simple  step  to  obtain 
R(m/N)3/2,  (H/H3/z,  and  5  apprQr>  ,  since 

r(m/h)3/'2  =  £  d/Vi3) 

(M/1D3/2  =  f£  (1/Vi2)]3^ 

and  £  approx.  is  Siven  by  (17)  or  by  a  six-term  formula  whose  coefficients 

are  those  of  (9).  Nothing  will  be  said  now  about  the  comparative  values 

of  the  three-  and  six-direction  values  of  £  .A  full  discussion  on 

approx. 

ot4  in  a  later  section  can  be  considered  to  be  quite  applicable  to  ®-aT)I)ro3: 
since  the  error  in  the  integral  cLft-  can  be  made  negligible  in 

comparison  with  that  of  ;  therefore  anything  said  about  the 

accuracy  of  the  Schiff  approximation  to  aC,  will  apply  to  ^8pprox  *  too. 


(a)  R(M/N)3/Z  x  io15  (cm/sec) 
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Determination  of 


j  (m/u) 


3/2 


iLa 


The  problem  remaining  is  to  determine  the  value  of  the  integral 

(  o  /  p 

l  (M/N)J  dA,  for  which  we  require  the  coefficients  of  the  cubic  equation 
in  v-  .  These  coefficients  are  obtained  in  terms  of  the  elastic  constants 
in  Appendix  A.  The  integral  becomes 


(20) 

J  (M/H)?/2  4A  =  (j>/eW)3/2 

=  (p/c/rt.)^2 

where 

p  =  3  +  2  (<r+  f  ) 

r  =  CT+  $  +  1 

s  =  (J2  (o-+3f  )  /4 

t  =  cr  («r+2S  ) 

and 

(21) 

CT  =  (cn  -  c12  -  2c^)  /c^ 

f  =  ^ci2  +  c44)  /°ii4  • 

f  and  g  are  functions  of  directions  and  are  given  by 

(22)  t  =  ij-x^z2 

=  cos*"  9  sin^  9  sin*1  2^ 

g  =  X2?2  +  y^z2  +  z2*2 

=  sin2  9  -  sin^  9  [l  -  (sin22*)/4] 

x,  7,  and  z  are  the  directions  cosines  of  the  wave-vector  k  whose  polar 

1 1 

co-ordinates  are  0,  , 

It  is  more  convenient  to  write  the  integrand  I  in  (20)  as 


(23) 


=  (I 


3/2 


ip  +  df  _ . 

3/2 

obtained  by  dividing  numerator  and  denominator  by  (  -  t  )J  . 
(24)  a  =-p/t 


Thus 


.  3+2  (<r+f ) 

'-c  (<r+  3  £ ) 


0  '  ■  0*1  *r  ■ 


. 


.  v 


. 


(  +  ■  - 
‘  ■  ;■  '■  -  ■: 

■ 

■-  •;  -  i 


.  ' 


-  *r 


\  >  “ 

. 


?v,  . 
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c  =  -  r/t  =  <T  +  £+  \ 

-<r (<r+  2$  ) 

d = -  s/t  =  -  ^£±4- 

For  all  cubic  metals  V  <  0  and  for  elastic  stability,  it  is  necessary 
that  S'  >  -2.  In  general  «r  decreases  with  increasing  anisotrox>y  and  is 
zero  for  isotropic  substances.  For  all  cubic  metals  considered  in  this  work, 
the  integrand  I  is  real  and  finite  over  the  whole  un.it  sphere. 

V/e  will  now  consider  three  methods  of  evaluating  (20),  each  of  which, 
through  reasons  of  accurancy  or  amount  of  labour  involved,  is  well-suited 
to  only  a  portion  of  the  <T  interval. 

Method  I.  Taylor  Expansion  About  f  and  g. 

V/e  expand  the  integrand  I  in  a  Taylor  series  about  f  =  4/ 10 5  and 
g  =  1/5,  which  are  their  average  values  over  the  unit  sphere. 

(25)  I  =  [I]  +  gi]  (f  -  4/105)  +  gl]  (g  -  1/5) 

+  M  Pf^|  (f  "  4/105) 2  +  2  j^lj  (f  -  4/105)  (g  -  1/5) 

{ -friz]  (s  “  ll5)2]  +  •  •  • 

where  [x]  indicates  V  taken  for  average  values  of  f  and  g. 

In  those  cases  where  the  crystal  is  nearly  isotropic  (  e.g.,  W  and  Al) , 
terms  in  f  may  be  regarded  as  negligible  in  comparison  with  those  in  g  only, 
since  d,  the  coefficient  of  f,  decreases  with  increasing  isotrapy. 

The  Taylor  series  is  integrated  term  by  term  over  the  surface  of 
the  unit  sphere. 

(26)  (1/4  tt)  j  I  cm  =  [1]  +  j^]  (1/4  it)  j  (f  -  4/105)  d-CL 

+  (1/4  tt)  |  (e  -  1/5)  dUl  +  .  .  . 
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_  ,  con5is4s 

Each  term  -  of 


(a)  an  integral  (l/4n)  J  (f  -  4/l05)p 


(g  -  l/5)q  drt.. 


and  Cb)  [xj  whose  value  depends  on  the  elastic  constants  of  the  material. 
The  integral  (l/4n)  j  (f  -  4/l05)p  (g  -  1/5)  q  cLa  is  tabulated  helm/ 
for  values  of  p  and  q_  up  to  5» 


TABLE  4, 


also  require  the  values  of  J  ,  taken  at  f  =  4/105  and 

General  expressions  for  (  d n  I)/  OfP  dgQL)  -up  to  n  =  5  (p  +  q_  =  n) 
found  in  Appendix  B.  These  are  applied  to  W,  Al,  Ee,  and  Pb;  the 
a,re  tabulated  below.  Following  the  notation  of  the  Appendix,  we 


We 

g  =  1/5 

v/ill  be 
results 
write 

i  =  (P/q) 

where  P  =  a  -  g 

Q,  =  c  +  df  - 


3/2 


nT> 


5  / 


, 

■'  c  _r'-  ■  .  .  -rx:?o  do.rh 

'  _  -  '  ‘  '  : 
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•  _  •  .  ' '  '  -  ■ :  v  ■  "  -■ 
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TABLE  5. 

(dnD/(afPdg^) 

(a)  Tungsten 

P  =  -249.  28 
Q  =  -108.43 
[l]=  (P/Q)3/2  =  3.4857 

The  other  derivatives  turn  out  to  he  negligible  in  comparison,  and  are 
therefore  neglected.  In  the  tables  below,  the  maximum  value  of  p  or 
is  the  same  as  the  number  of  terms  taken  in  the  Taylor  expansion. 

(b)  Aluminum 


p  = 

4.3989 

q  =  1. 8474 

d.  =  0.11843 

X 

0 

1  2 

3 

0 

3.6744 

-0.35331  0.056626 

-0.012705 

1 

1.7304 

-0.35765  0.087969 

2 

2. 1451 

-0.59343 

3 

3.8817 

(c)  Iron 


P  =  1.1313  Q,  =  O.35356  d  =  0.49238 


•  . 

. 

. 

!  ■  .  ■J-  ;• '  *  ’  ••• 

.  :  '  •;  ;  •>  - :  -  ...  ••  Ho  ' : .  ss:.  '■  \  -• 
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(&)  Lead 


P  =  0.62397  Q,  =  0.18266  d  =  0.60674 


All  that  now  remains  to  obtain  (i/4tt)  j  XcLa.  is  to  combine  the 
quantities  of  Tables  4  and  5  according  to  (26) .  For  purposes  of 
indicating  the  degree  of  convergence,  the  individual  terms  of  the 
series  are  tabulated  below. 


SABLE  6. 


Metal 

Term  0 

Term  1 

Term  2  Term  3 

Term  4 

Term  5 

Total 

¥ 

3.4857 

0 

0 

3.4857 

A1 

3.6744 

0 

0.0072  -0.0002 

3.6814 

Fe 

5.7237 

0 

0.2676  -0.0569 

0.0355 

-0.0142 

5.9557 

Pb 

6.3137 

0 

1.0446  -0.4288 

0.4973 

-0.3804 

7.0464 

The  number  of  terms  required  depends  on  the  anisotropy  of  the 
crystal.  Since  the  labour  in  calculating  the  terms  is  quite  considerable, 
no  attempts  was  made  to  go  beyond  five  terms.  It  was  found  that  five 
terms  were  sufficient  to  calculate  the  integral  to  better  than  1$  (in 
comparison  with  accurate  results  using  elliptic  integrals)  \fhen  <r  was 


greater  than  -1 


-  ■  ■» -j*  ■■ 


?  •  ■  ■  •  .  ;  !•  ; 
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Method  II.  Taylor  Expansion  about  sin^2^=  1/2 


Consider  the  relative  change  in  direction  of  a  unit  vector  due  to 
changes  in  ©  and  if.  It  is  obvious  that  over  a  large  portion  of  the  unit 
sphere,  the  direction  does  not  change  as  much  for  a  change  in  only  as 
it  does  for  one  in  ©  only  (this  is  especially  evident  for  small  ©)• 
Consequently,  the  function  of  direction  in  the  integrand  will  not  vary 
in  general  as  much  with  ij)  as  with  ©,  so  that  an  expansion  of  the 
integrand  with  respect  to  cj)  and  a  numerical  integration  with  respect  to 
©  should  he  an  improvement  over  Method  I,  even  when  the  number  of  terms  in 
the  expansion  is  small.  Since  the  integrand  is  a  function  of  through 

2  p 

sin  2  if,  we  expand  about  sin^2ij?=  1/2.  The  integrand,  in  terms  of  ©  and 
f,  is 


(28)  I 


a  -  sin2©  +  sin1©  (l  -(sin22tp)  /4) 


3/2 


.  c  +  d  cos2©  sin^©  sin22(f-  sin2©  +  sin^©  (l  -  (sin22<f)  jk\ 


where 


Thus, 


(29)  I  =  [i]  + 


with  [x]  indicating  that  x  is  calculated  at  sin22i{?=  l/2.  Integrating 
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the  separate  terms  with  respect  to  and  simplifying,  we  obtain 

(30)  d/at)  f  "i«.  +  2  M  -wzj,(xy. jwz,  r.g&  +  . 

''O  1  \Y  +  Z/2 1  64  (W  +  X/ 2)1/,2(Y  +  Z/2)7'2 

The  G  integration  is  performed  numerically,  taking  intervals  of  cos  G  =  0 

The  use  of  only  the  two  terms  in  (30)  is  sufficient  to  give  results 

accurate  to  within  1 fa  forcr>  -1.5.  Except  for  c r  values  near  0,  a 

region  where  Method  I  is  best,  this  method  appears  to  be  the  quickest 

of  the  three  discussed  here. 

We  apply  this  method  to  the  calculation  of  j" Idl2 for  Ee,  Ag,  Cu, 
and  Pb.  The  quantities  W,  X,  Y,  and  Z,  calculated  at  intervals  of 
cos  G  =  0.1,  appear  in  Appendix  C.  In  Table  7,  we  combine  them 
according  to  (3$)*»  final 3.y,  we  obtain  (l/4rr)  j*  Idilby  integrating 
numerically  by  Simpson5  3  rule. 


TABLE  7. 


(a)  Iron 


cos  G 

Term  1 

Term  2 

(1/2 n) 

J 

Air 

IdcP 

0  1 

(l/4rr) 

IdQ. 

0.0 

3.0507 

0.2711 

5. 3218 

0.1 

5. 0971 

0. 2568 

5; 3539 

0.2 

5. 2422 

0.2159 

5.4581 

0.3 

5.5013 

0.1556 

5. 6569 

0.4 

5.8891 

0.0876 

5.9767 

0.5 

6.3945 

0.0298 

6.4243 

0.6 

6.9094 

0.0007 

6.9051 

0.7 

7.1073 

0.0004 

7.1077 

0.8 

6.5898 

0.0029 

6.5927 

0.9 

5.35°^ 

0.0004 

5. 3508 

1.0 

3. 9277 

0.0000 

3. 9277 

5. 9563 
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(t>)  Silver 


cos  0 

Term  1 

Term  2 

(l/2n) 

(i/4tt) 

ua 

0.0 

4. 9614 

0.4535 

5.4149 

0.1 

5.0158 

0.4312 

5.4470 

0.2 

5.1889 

0.3672 

5.5561 

0.3 

5.5099 

0. 2699 

5.7798 

0.4 

6.0230 

0.1549 

6.1779 

0.5 

6.7491 

0.0501 

6.7992 

0.6 

7.5560 

-0.0025 

7.5535 

o.7 

7.9139 

0.0064 

7.9203 

0.8 

7.1000 

0.0103 

7.1103 

0.9 

5.3316 

0.0010 

5.3326 

1.0 

3.6355 

0.0000 

3.6355 

6. 2321 

(c)  Copper 


cos  Q 

Term  1 

Term  2 

(l/2n) 

J 

rin 

Mf 

(l/4rr) 

Idil 

0.0 

5.3253 

0.5639 

5.8892 

0.1 

5.3841 

0.5265 

5.9106 

0.2 

5.5733 

0.4435 

6.0168 

0.3 

5. 9281 

0.3187 

6. 2468 

0.4 

6.5039 

0.1742 

6. 6781 

0.5 

7.3353 

0.0478 

7.3831 

0.6 

8.2897 

-O.OO55 

8. 2842 

0.7 

8.7437 

0.0189 

8.7617 

0.8 

7.8037 

0.0186 

7.8233 

0.9 

5.7427 

0.0016 

5.7443 

1.0 

3.8171 

0.0000 

3.8171 

6.7832 

, 


(o) 

“ 


-SO- 


ld)  Lead 


cos  © 

Term  1 

Term  2 

(I/3n) 

"a.ir 

(1/4ti) 

Id  11 

0.0 

5.1445 

0.8491 

5.9936 

0.1 

5.2101 

0.8094 

6.0195 

0.2 

5.4247 

0.6947 

6.1194 

0.3 

5.8468 

0.5172 

6.3640 

0.4 

6.5834 

0.2972 

6.8806 

0.5 

7.7633 

0,0816 

7.8449 

0.6 

9 . 2956 

-0.0123 

9.2833 

0.7 

10.1004 

0.0697 

10.1701 

0.8 

8.5475 

0.0531 

8. 6006 

0.9 

5.6417 

0.0032 

5.6449 

1.0 

3.4691 

0.0000 

3.4691 

7.1801 

Method  III,  Application  of  j£Lliptio  Integrals  to  <|)  Integration 

For  cr>  -i.5,  we  must,  however,  resort  to  the  more  accurate,  as 
well  as  the  most  labourious  method,  wherein  the  cj)  integration  is 
performed  in  terms  of  elliptic  functions,  followed  by  numerical 
integration  with  respect  to  ©.  Let  us  write  (28)  as 


(31) 


1  r  (W/Y) 


3/2/1 


-  n|sin22<p  \5/2 


1  -  n^sin  2cf 


nx  -  -  X/W 


4 

sin  © 


4(a  -  sin^©  +  sin4©) 


l2  -  -z/Y  r  ^in4Q  .  (1/4  -  d  cos^©) 


c  -  sin2©  +  sin4© 


where 


. 


X. 


* 


. 


.... 

. 
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ITow  consider  the  integral 
(3E) 


(W/Y)3/2 

Cl 

a  - 

nfsin22(pf/2 

dcp 

V 

,i  - 

“  "*J  P  "4  "  j 

■  ngSin  2cp  / 

1 

=  t(W/Y) 

3/2  | 

pi 

f  1  - 

nfsinV\ 

3/2 

| 

dr¬ 

u  - 

nfsin V  ) 

=  4(W/Y) 

3/2 

r  it/* 

n  - 

n2sin*V  ^ 

O  Q 

3/1 

| 

ay 

o  1 

\  1  - 

n*sin“v  y 
&  • 

z  4(W/Y) 

3/2  • 

I. 

As  the  integration  of  J  is  dependent  on  the  absolute  and  relative 

o  g 

magnitudes  of  n£  and  n2,  it  must  he  independently  performed  for  the 
three  regions 

X.  n3  <  0 

II.  1  >  nf  >  n|  >  0 

III.  1  >  n§  >  nf  >  0 

The  integration  is  carried  out  in  Appendix  D,  and  only  the  results  will 

1  £3 

be  quoted  here*  The  notation  is  that  of  Byrd  and  Friedman  • 

(33a)  I. 

r  2n 


4  <  0 


Jo 


idcf  -  (W/Y)3/3  {4k'S  +  6nfk,3(E  -  k,2K)/k2  -  6nfk'3Kj 


with  k3  =  1  -  k'  =  |n^(l  +  In,,  |) 


(33b) 


II. 


1  >  nf  >  ni  >  0 


I4<p  =(WA) 3/2  j  4g3  -  *4  rgE  -  1  /Yj-ijr, 

3  L  Ho  L  2 


;,k) 


with  k3  =  1  -  k>2  =  ?1, ;  g3  =  1/(1  -  nf)  ;  iir  =  sin"1 

-r  2  T  _ 


1  -  n 


A  0(^>k)  =  (2/tt)  |EP(ijr,k’)  f  KEOjr.k*)  -  KS’(ljr,k,)| 
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III.  1  >  n|  >  >  0 


Tables  of  the  functions  K,  E,  andA0»  well- suited  to  the  evaluation 
of  these  integrals,  are  available  in  an  appendix  in  Ref,  18,  Interpolation 
is,  however,  almost  always  necessary. 

Method  III  is  applied  to  the  calculation  of  J Id!) for  the  fairly 
anisotropic  metals  Ag  and  Cu  and  the  very  anisotropic  metals  Pb,  K,  and 
Ha.  Most  of  the  numerical  work  is  tabulated  in  Appendix  E;  however, 


[  Id  Cp  and  (l/4ir)  fldfihave 

^  fJkTr 

J  IdCp-  is 


the  quantities  of  interest,  namely  (1/2tt) 

been  gathered  in  Table  8,  In  addition,  the  quantity  (1/2tt) 
plotted  against  cos  Q  for  Ag,  Pb,  and  Ha  in  Pig,  1.  The  effect  of 
anisdtropy  is  very  evidait.  We  note  especially  the  increasing 
contribution  of  the  peak  as  we  go  from  Ag  to  He,. 
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TABLE  8. 


cos  Q 

Ag 

Cu 

Pb 

K 

Na 

0.0 

5.4945 

5.9917 

6.3038 

8.7403 

11.1353 

0.1 

5.5174 

6.0089 

6.2939 

8.6319 

10.7949 

0.2 

5.6051 

6.0835 

6.3094 

8.4244 

10.1100 

0.3- 

5.8052 

6.2793 

6.4598 

8.3926 

9 • 6889 

0.4 

6.1855 

6.6874 

6.9102 

8.8563 

10.0611 

0.5 

6.7997 

7.3842 

7.8460 

10.1808 

11.8388 

0.6 

7.5514 

8.2835 

9.2810 

12.5325 

15.7185 

0.7 

7.9113 

8.7519 

11.6428 

14.4444 

19.6127 

0.8 

7.1086 

7.8204 

8.5971 

11.9334 

15.0013 

0.9 

5.3323 

5.7441 

5.6449 

7.2502 

8.0127 

1.0 

3.6351 

3.8169 

3.4693 

4.1084 

4.1875 

ir/MA- 

6.2498 

6.8077 

7.2534 

9.7313 

11.897 

Probable  Error  in  Method  III. 

Since  the  criterion  of  accuracy  for  Methods  I  and  II  is  Method  III, 
it  is  necessary  that  we  have  sajie  idea  of  the  magnitude  of  possible 
errors  in  III.  The  sources  of  error  in  this  method  are  two  : 

(1)  the  numerical  Q  integration,  and 

(2)  the  use  of  a  two- tern  binomial  expansion  for  (1  -  n?sin  V j 
in  region  I.  (see  Appendix  D) . 

It  turns  out  for  the  metals  considered  that  n^  is  quite  small  in 
this  region,  never  exceeding  0.2  for  the  alkali  metals  and  is  less  than 
metals,  (see  Appendix  E) .  When  average  values  of  n^  and 
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Fig.  1. 


Zrr 

I  dy  by  Method  III  for  Ag,  Pb,  and  Na. 
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2  °  2 

sin  Y  are  introduced,  one  finds  that  1  -  (3/2)n£sin  Y  differs  from 

(1  -  n^sin2)Q^/3  by  about  1/10  of  1%  for  Na.  The  actual  error,  of 
course,  is  even  less,  since  Region  I  represents  about  1/3  of  the  whole 
range  of  integration. 

The  error  introduced  by  the  numerical  integration,  though  greater 
than  that  due  to  the  previous  source,  is  still  quite  small.  A  rough 
cheek  for  sodium  using  integration  intervals  of  0.05 gave  a  result  that 
differed  by  about  1/4  of  lf5  from  that  obtained  using  intervals  of  0.1. 
It  would  appear  safe  to  conclude  that  the  accuracy  of  the  method  is 
about  1/3  of  1 fa  for  highly  anisotropic  metals  and  improves  with 
increasing  isotropy.  It  is  certainly  valid  to  use  it  as  a  measure  of 
accuracy  of  Methods  I  and  II. 

( 1/ 4tt)  j  Idfl  by  the  Three  Methods 

To  facilitate  comparison,  we  collect  the  results  for  (1/4tt)  J  Idji 
obtained  from  the  three  methods  into  a  table. 


TABLE  9. 


Metal 

or 

I 

II 

III 

W 

0.0066225 

3.4857 

A1 

-0.35714 

3.6814 

Fe 

-1.1724 

5.9557 

5.9563 

Ag 

-1.3066 

6.2321 

6.2498 

Gu 

-1.3750 

6.7832 

6.8077 

Pb 

-1.5000 

7.0464 

7.1801 

7.2534 

K 

-1.6538 

9.7313 

Na 

-1.7586 

11.897 
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From  the  table,  and  keeping  in  mind  previous  comments  on  the 
relative  amount  of  labour  involved  in  the  application  of  these  methods, 
it  would  appear  that  Method  II,  despite  the  fact  that  it  is  less  accurate 
than  III,  is  quite  suitable  for  most  purposes  for  O' values  in  the  range 
-0.5  to  -1.5.  Above  CT  -  -0.5,  we  find  Method  I  to  be  sufficiently 
accurate  as  well  as  quickly  applied.  Below  O’  s  -1.5,  it  is  necessary 
to  use  Method  III  if  the  Debye  9  accurate  to  3$  is  desired. 

Comparison  of  the  Bhatia  -  Tauber  and  Schiff  Approximations  toflC-. 

’The  data  in  Table  3(a)  can  be  applied  with  little  extra  effort  to 
the  calculation  of  by  the  B  -  T  method,  since  the  R(M/U}3^2  need  only 
be  added  after  multiplication  by  the  appropriate  coefficients.  We  will 
denote  the  various  approximations  as  : 

A.l  -  three- term  Bhatia  -  Tauber  approximation,  as  given  by  (8)  ; 

A. 2  -  six-term  Bhatia  -  Tauber  approximation  as  given  by  (9)  ; 

A. 3  -  Schiff  approximation  using  three- term  R.  di  is  given  by 
«  (l/2Tr2)l“(l/4ir)  JldJX  ; 

A. 4  -  Schiff  approximation  using  six-term  R. 

The  results  obtained  from  the  four  approximations  are  compiled 
in  Table  10  on  the  next  page. 


;  c 
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TABLE  10  .  c<i  x  1018 


Metal 

B  -  T  Method 

Sehiff  Approximation 

A.l 

A. 2 

A.3 

A.  4 

W 

4.8880 

4.8880 

Method  I  i 

A1 

3.0652 

3.056 

I 

Fe 

3. 6656 

3.470 

II 

Ag 

27.736 

24.452 

25.257 

24.614 

III 

Gu 

10.692 

9.2483 

9.6093 

9.3204 

III 

Pb 

231.44 

181.56 

192.  78 

184.26 

III 

K  . 

91.900 

64.789 

70.718 

66.083 

III 

Fa 

51.870 

30.446- 

34.292 

31.286 

III 

From  theoretical  considerations,  one  would  expect  the  Sehiff 
approximations  (A, 5  and  A. 4}  to  lead  to  better  results  than  those  of 
Bhatia  -  Tauber  (A.l  and  A* 2),  and  this  seems  to  be  confirmed  by- 
computations  for  actual  eases.  It  can  be  seen  that  calculations  using 
A. 4  do  not  differ  substantially  from  those  of  A. 3,  indicating  that 
A. 3  is  already  very  good.  Therefore,  one  expects  that  A. 4  leads  to 
very  nearly  the  correct  results.  This  is  not  the  behaviour,  however, 
of  the  Bhatia  -  Tauber  approximation,  where  the  great  improvement  of 
A. 2  over  A.l  suggests  that  A. 2  itself  could  be  improved,  perhaps 
considerably,  by  using  more  than  six  terns  in  the  B  -  T  formula. 

A  comparison  of  as  calculated  by  A. 2  and  A. 4  does  indicate  that 
A.  2  is  indeed  quite  good  ;  in  fact,  it  is  within  3$  (1>j  in  9  )  of 
A. 4  for  all  or  values  down  to  -1.75  and  this  region  contains  perhaps 
nearly  all  of  the  cubic  metals.  A  point  in  favour  of  A, 2  is  that  it 
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involves  only  one- third  to  one-half  of  the  calculations  required  in  A* 4, 

The  accuracy  of  A. 4  depends  strictly  on  that  of  R,  since  J  IcLQ.  can 
be  obtained  almost  as  exactly  as  required.  Without  calculating  the 
influence  of  higher  Kubic  Earmonics  on  B,  it  is  not  possible  to  state 
definitely  how  good  R  is,  but  considering  that  A. 4  represents  an 
improvement  over  A. 3  of  10$  for  the  worst  case  (Na),  it  is  probable 
that  A. 4  leads  to  results  for  «r,  that  are  not  out  by  more  than  a  couple 
of  percent.  Thus  0^  by  A. 4  is  probably  better  than  1$  in  the  case  of 
Na,  with  even  better  results,  of  course,  for  more  isotropic  metals. 

Results  for  9^ 

Finally,  using  the  beSt  values  of  Jldikand  R,  the  following  results 
for  have  been  obtained. 


TABL33  11 


Het  al  ; 

R 

(l/4w)  flciit 

G2°E 

This  Work 

Blackman^9 

Q,uimby° 

W 

0.62015 

3.4857 

371 

Al 

0.63512 

3.6814 

429 

!e 

0.64993 

5.9563 

462 

461 

Ag 

0.66083 

6.2498 

212 

212 

211.3 

Gu 

0.66015 

6.8077 

332 

Pb 

0.68785 

7.2534 

89.5 

K 

0.89251 

9.7313 

93.4 

Na 

0,71802 

11,897 

150 
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Comparisons  between  these  results  and  any  experimental,  values 
that  might  be  available  are,  with  two  exceptions,  not  valid,  since 
the  elastic  constants  data  used  here  were  obtained  at  room  temperature. 

In  addition,  the  use  of  different  constants  makes  it  useless  to  quote 
the  results  of  other  theoretical  calculations  available  in  the  literature. 
However,  in  the  case  of*  Fe,  Ag,  and  Ha,  comparisons  with  values  found  in 
Blackman‘S  are  possible,  since  the  elastic  constants  used  there  are  the 
same  (or  very  nearly  so),  as  those  c£  this  work*  In  addition,  9^  -Lor  silver 
calculated  on  the  basis  of  the  Q,uimby  -  Sutton  approx  mat!  on  is 
available  for  comparison. 

Accuracy  of  Calculations 

On  the  basis  of  the  approximations  used,  the  accuracy  of  the  results 
for  9^  is  estimated  to  be  better  than  It  is  not  likely  that  the 

numerical  work  would  introduce  any  significant  errors,  since  the 
calculations  were  almost  always  carried  out  using  five  figures.  It  is 
to  be  understood,  however,  that  only  three  of  these,  at  most,  are 
significant,  as  the  experimental  data  used  are  not  usually  known  any 
better  than  that.  It  would  seem  safe  to  conclude  then,  that  the  Schiff 
approximation  A. 4  is  at  least  as  good  as  any  of  the  data  that  one  would 
generally  have  to  use. 
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PART  III.  DEBYE  CHARACTERIST IC  TEMPERATURE  OP  SILTER 
General  Expression  for  the  Specif ic  Heat 

The  specific  heat,  cv,  per  atom  at  constant  volume  is  given  by 
(I. 10)  as 

(1)  H0Vcv  =  VE(ti/KT)2  jff(ra)ai2eIlIl,,/KT(elm/KT  _  l)“2do, 
where 

(2)  ygM&v  =  S-  a4  >  k2(®.,G,<<>)  d_n_ 

\2n)~>  J  fz 7  1  du^ 

Approximations  to  G-(o)) 

Por  the  determination  of  G(co)  we  require  the  complete  dispersion 
relations  "between  o>  and  k.  It  is  clear  that  what  has  been  said  in  Part  II 
with  regard  to  the  complexity  of  these  relations  in  the  long  wavelength 
limit  and  to  the  difficulty  of  carrying  out  the  necessary  integration 
obtains  here  fortiore.  We  thus  again  resort  to  approximate  methods. 

The  two  methods  used  here  are  (l)  the  three-term  Houston  approximation, 
and  (2)  and  alternative  three-term  formula  whose  coefficients  are 
determined  by  the  requirements  that 

(a)  Gr(u))  give  the  correct  number  of  modes,  and 

(b)  at  T  =  0°K  agree  with  the  experimental  value. 

Houston  Approximation’ 

We  have  already  encountered  the  Houston  approximation  in  Part  II, 
so  that  it  is  sufficient  here  to  merely  write  down  the  expression  for 
G-(uj)  obtained  from  it. 

(3)  &(o>)  =  llO^dk/aro^  +  l6(k2dk/ dap-g  +  9(k2dk/to.)c| 
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The  factor  F  is  introduced,  to  ensure  that 

(4)  JgWcLio  =  3Uq. 

Determination  of  F 

If  we  denote  the  maximum  value  of  k  along  the  directions  A,  B,  and 
C  byk^(in),  k^(m),  and  kg (m),  respectively,  then 

(5)  Jg-(u>)6/i)  =  |l0k£(m)  +  l6kg(m)  4-  9kg  (m)j 

=  3^ 

The  k-vectors  are  confined  to  the  first  Brillouin  zone,  since  solutions 
for  k  outside  this  zone  give  no  new  modes  of  vibration.  Therefore  the 
k(m)  are  to  he  determined  by  finding  the  point  of  intersection  of  the 
k-vector  with  the  boundary  of  the  first  Brillouin  zone.  For  fee  crystals, 
the  first  Brillouin  zone  is  marked  out  by  the  (200)  and  (ill)  planed  and 
the  k(m)  are20 

(6)  k^(m)  =  2(*rr/a) 

kg  (m)  =  3.(Tr/ a) 

kg(m) 

Squat ions  (5)  and  (6)  give  for  F: 

(7)  F=  0.95664. 

Physical  Basis  of  F 

Physically,  F  arises  from  the  deviations  of  the  Brillouin  zone  from 
spherical  symmetry.  For  small  k,  whose  contributions  are  significant  at 
low  temperatures,  the  shape  of  the  Brillouin  zone  is  clearly  of  no 
importance.  Thus,  no  F  factor  is  required  in  this  case,.  As  the 
temperature  rises,  ever  larger  k  values  contribute  to  the  specific 
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heat,  until  k  values  of  the  magnitude  of  k(m)  are  significant.  However, 

these  k(m)  are  not  of  equal  length  in  every  direction,  and  the  resulting 

asymmetry  plus  the  use  of  a  finite  number  of  directions  in  Houston* s 

approximation  result  in  a  necessity  for  an  adjusting  factor  F.  The  F 

that  is  used  in  the  calculations  is  an  average  value,  since  no  factor  is 

required  for  small  k,  and  some  factor  larger  than  F  is  required  in  the  region 

beyond  the  Brillouin  surface  closest  to  the  origin.  Thus,  the  introduction 

of  this  factor  into  G-(u))  cannot  in  general  lead  to  the  right  results 

for  cy.  at  low  temperatures.  However,  in  the  case  of  silver,  the  0-p  at 

0°K  obtained  by  this  method  (i.  e.  the  B  -  T  formula)  is  in  agreement 

21 

with  the  latest  experimental  value.  This  is  certainly  due  to  a  number 
of  compensating  factors:  (l)  the  use  of  a  three- term  approximate 
formula,  (2)  the  use  of  approximate  values  for  the  elastic  constants 
at  0°H,  and  (3)  the  possibility  that  the  experimental  results  are  not  exact. 

Adjusted  Constants  Method 

In  the  second  method  of  calculating  the  frequency  spectrum,  the 
Houston  coefficients  10,  1 6,  and  9  are  replaced  by  16,  a,  and  b,  and 
a  and  b  are  then  determined  by  requiring  that 


(8) 


where  otq  is  the  value  obtained  by  A. 4  (Part  II).  These  two  conditions 
automatically  ensure  that  0^  is  correct  at  the  two  extremes  of  low  and 
high  temperatures.  Good  results  at  the  high  temperature  end  are  assured 


because  | &(u))E(x)du)  is  nearly  independent  of  the  shape  of  the  G(u))  curve, 
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since  E(x)  changes  very  slowly  with  u)  when  T  is  large  (x  =  tiu)/KT).  What 
we  nearly  have  is  a  constant  x  Jg-(o))&u>  and  j'Cr(u>)^has  already  been  fixed 
hy  the  first  condition*  At  the  low  temperature  end,  agreement  is  obtained 
through  This  method  will  hereafter  he  referred  to  as  the  ’’Adjusted 

Constants  Method11  (AC  Method)* 

Determination  of  a  and  b 

Because  it  is  the  most  accurate  value  of  oC^  that  we  have,  we  will  use 
the  one  obtained  from  A. 4  (Part  II)  for  determining  a  and  b.  However, 
before  inserting  it  in  (8«4>  must  first  convert  it  to  the  0°K  value, 
since  it  was  calculated  using  room  temperature  elastic  constants,  the 
only  temperature  at  which  the  elastic  constants  for  silver  have  been 
measured.  In  a  later  section  our  approximate  extrapolation  of  the  room 
temperature  values  of  the  elastic  constants  to  0°K  is  equivalent  to 
multiplying  by  (I.I83 )~3/2  =  0.7772  to  convert  it  to  the  0°K  value 
to  give 

(9)  =  19.15°  2:  10"1®  (cm/ sec)~3  (0°K) 

Then,  from  (8)  we  obtain  (using  IT  =  5*9  x  1022/Cc), 

(10)  a  =  7.2083 
b  =  7.1006 

and  the  approximate  expression  for  the  frequency  spectrum  becomes 


(11)  G(®)  = — l—r-JB—  y  jioCk^dk/auOi  +  7.2083(fc2ak/aa)i)B 

(2n)J  24.309  7T,  '• 

+  7.1006(k2dk/d!ui)c| 

In  comparing  (ll)  with  the  Houston  approximation  (3),  we  note  that 
the  contribution  from  the  B  direction  has  been  considerably  decreased 
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while  the  contributions  from  the  A  and  C  (especially  A)  directions  have 
been  increased. 

Dispersion  Relations 

In  order  to  obtain  Gr(uu)  by  either  of  the  approximations  (3)  or  (ll) 
we  require  the  solutions  of  the  secular  equation  (1.28)  along  the  three 
principal  directions  of  the  crystal,  A,  B,  and  C.  They  are: 

(12a)  Direction  A  (100) 

=  (8c  11  / a2p  )sin2  (ak/  2/2) 

®tl  =  ^8c44  /a2p  )sin2  (ak/ 2^2) 

~  ^sin2  (sk/2j2) 

(l2b)  Direction  B  (110) 

=  (8/a,2 p  )  sin2  (ak/4)  j  2e1]L  -  €  -  (20-q  -  -  *  )sin2  (ak/4)j 

u.f  =  (8/a2p  )sin2  (ak/4)  [  20^  -  (2044  -  cn)  sin2(ak/4)} 

u)|2  =  (8/a 2p  )sin2  (ak/4)  {  2044  +  (•  -  (044  +  fe  )sin2  (ak/4)j 

(l2c)  Direction  C  (ill) 

ot>2  =  (2/a' p  )  -  2  *)sine(ak/ /~6) 

w^l  =  (2/a2p  )  (3e^  +  €  )sina(ak j J6) 

a)2  =  {z/aFp  )  (3c  +  €■  )sin  e  (ak//6) 

t2  r  44 

Temperature  Dependence  of  Gr(u>) 

Because  of  the  dependence  of  the  elastic  constants  on  temperature,  it 
is  that  Gr(a>)  will  vary  with  temperature.  Thus,  strictly  speaking, 

in  calculating  cy  at  any  temperature,  the  Gr(o))  appropriate  to  that  temperature 
should  be  used.  It  can  be  appreciated  that  the  labour  involved  in  carrying 
out  a  complete  calculation  of  the  Debye  temperature  over  the  whole  temperature 
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range,  taking  into  recount  the  detailed  variation  of  the  elastic  constants, 
is  quite  substantial.  However,  for  the  case  of  silver,  some  simplifying 
assumptions  can,  of  necessity,  be  made,  although  it  is  known  that  they  are 
not  strictly  true.  They  are: 

(1)  that  the  elastic  constants  vary  linearly  with  temperature,  and 

(2)  that  the  t  emperature  coefficient  of  the  constants  is  the  same 
for  all  of  them. 


Adjustment  of  G-(ou)  for  Temperature  Variation 

Assuming,  then,  that  the  three  elastic  constants  vary  with  temperature 
in  precisely  the  same  way,  it  is  clear  from  (12)  that  all  the  dispersion 
relations  can  be  written  in  the  form 
03)  u)  =  °441/2  f  (k) 


and  that  the  temperature  dependence  of  uj  is  contained  solely  in  the 

i/2 

(f°r  convenience,  we  will  leave  out  the  subscript  44  in  what 

>  1/2 

Let  x  =  w/c  •  then  we  can  write 

at)  kz 

-  c-1/2  k2  f§ 


factor 
follows) . 


2  dk 

where  k  -=•  is  some  function  of  x.  Thus 
dx 

(15)  &(x)  =s  c  ^  g(x) 

where  g(x)  involves  a  summation  of  the  k^  ~~~  terms.  Denoting  at 

o  * 

0  K  and  at  any  other  temperature  T  by  c0  and  c^,  respectively,  we 

for  a  given  x 

(16)  gW  =  C^Z  Q.0(x)  =  «*T1/2STW 

Hence 

(17)  GjU)  =  (co/ct)1/2Go(x) 


have 
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and 

(18)  Gt(<o)  =  (co/cT)l/2  Gq  [^/e^,)172^] 

Thus,  once  the  frequency  spectrum  is  known  at  0°K,  it  can  easily  he 
determined  at  any  other  temperature. 

Choice  of  Elastic  Constants  for  0°K 

As  previously  mentioned,  measured  values  of  the  elastic  constants  for 
silver  are  known  at  present  only  for  room  temperature.  Thus,  to  determine 
the  C°K  values, it  is  necessary  to  carry  out  some  reasonable  extrapolation 
procedure.  Since  the  results  obtained  here  are  to  be  compared  to  those  of 
Bhatia  and  Horton^  ,  their  procedure  is  employed  here,  too.  The  room 
temperature  values 

(19)  cn  sr  12,0  ) 

°12  =  )  x  10U  dynes/ cm2  (T  =  .29l°K) 

) 

were  converted  to  0°K  values  by  multiplying  the  constants  by  the  factor 
1,183,  giving 

(20)  c,,  =  14.19  ) 

) 

a.,  =  10.62  )  x  1011  dynes /cm2  (5  =  0°K) 

12  } 

c44  -  5»17  ) 

If  we  calculate  <*,  by  A. 4  (Part  JI  )  using  the  elastic  constants  as 
given  by  (20) ,  we  obtain  (see  eq>  9) 

=  19.150  X  10"1S  (om/seo)-3  • 

How  this  value  of  eC|  leads  to  a  result  for  at  0°K  that  is  about  6° 
higher  than  the  latest  experimental  value^1  .  Since  A. 4  is  known  to  be 
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very  nearly  exact,  we  conclude  that  our  extrapolation  procedure  for  the 
elastic  constants  is  probably  not  good  enough.  In  an  approximate  attempt 
to  correct  this  discrepancy  between  the  theoretical  and  experimental  values 
of  ©D,  we  will  use  an  extrapolation  factor  1.123  instead  of  1.183.  Although 
it  may  first  appear  that  this  would  require  a  recalculation  of  the 
coefficients  in  our  approximation  (11)  for  G(w) ,  it  can  clearly  be 
incorporated  in  (18).  The  elastic  constants  determined  on  the  basis  of  the 
extrapolation  factor  1.123  are 


( 21)  c44  -  13.48 

c^ g  =  10.07 


x  1011  dynes/cm* 2.  (T  =  0  °K) 


c44  *  4.906  J 

Other  Constants  Used  in  Calculations 

The  other  constants  needed  to  calculate  Gr(co)  are  N0,  the  number  of 
atoms  per  unit  volume,  and  p,  the  density  of  the  substance.  These  are  taken 
to  be,  as  in  Ref.  2, 

(22)  N0  -  /2/a3  =  5.9  x  1022/cc  ;  £-10.5  gm/ce. 

and  p  also  vary  with  temperature,  but  this  variation  is  small  and  is  neglected. 


Procedure  Used  in  the  Calculation  of  G(  ) 

Ih  calculating  G(oj)  ,  rather  than  solving  for  k  and  dk/dco  f or  chosen 

values  of  co,  which,  at  least  in  the  case  of  Direction  B  would  be  very  difficult, 

the  reverse  procedure  is  used,  whereby  the  and  dw/dk  are  calculated  for 

2  2  l/2 

suitably  spaced  values  of  k.  A  k  dk/doo  vs.  (a  p)  curve  is  then  plotted 

for  each  solution  of  the  secular  equation  (seven  different  solutions  in  all 

2  l/2 

among  the  three  directions).  Finally, the  ordinates  corresponding  to  (a  ,g)  ‘ 
values  suitably  chosen  for  purposes  of  numerical  integration 
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are  combined  according  to  the  Hons ton  (3)  or  AC  (ll)  formula  to  give 
Cr(w).  For  small  values  of  co,  it  is  desirable  instead  to  use  the  analytic 
expression  for  &(oj)  ,  namely 
(23)  G(cd)  =  <aCjU)2 

=  Co)2a2 

since  this  region  contributes  significantly  to  cy  at  low  temperatures  and 

the  graphical  results  may  be  somewhat  unreliable  here.  Using  0°K  elastic 

constants  (20)  we  find  that  for  silver 

(24a)  &(u>)  =  2.4722  x  10“3  £  aV 

on  the  basis  of  the  Bhatia-Tauber  formula  including  F  and 

(241))  G(<d)  =  2.1938  x  10“ 3  p  a2®2 

on  the  basis  of  A.  4  (Part  II ). 

These  formulas  were  found  to  be  valid  to  lj6  for  ©(a2^)1^2  up  to  about  2  x  1CP* 
Treatment  of  Singularities  in  Gr(u>) 

Because  du)/<3 Ik  becomes  zero  when  k  =  k(m)  (except  for  the  two  transverse 
solutions  in  Direction  B) ,  the  frequency  distribution  function  as  determined 
by  either  of  the  two  methods  has  several  infinities.  However,  the  area 
under  the  curve  at  these  singularities  is  always  finite,  and  for  purposes 
of  numerical  integration,  the  singularities  can  be  replaced  by  finite  peaks 
of  equal  area. 

Comparison  of  G-(u>)  Obtained  from  the  Two  Approximations 

In  Fig,  2,  we  compare  the  &((»)  calculated  from  (3)  and  (ll).  The 
solid  curve  represents  (x(u>)  as  obtained  from  Houston* s  approximation  (3), 
while  the  dotted  curve  represents  that  obtained  from  the  AO  approximation  (ll). 
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We  draw  attention,  to  the  changes  in  the  relative  heights  -,:f 
of  the  peaks.  As  expected,  those  peaks  due  to  Direction 


the  contributions  to  G(a>)  from  Directions  A  and  C  has  an  appreciable  effect 
on  in  the  temperature  region  above  T  =  1 5°K,  but  this  will  be  discussed 
more  fully  in  a  succeeding  sections 

Debye  Temperature  of  Silver 

The  specific  heat  per  atom  at  constant  volume  can  now  be  calculated 
at  several  temperatures  from  (l)  using  numerical  integration.  Instead  of 
comparing  experimental  and  calculated  values  of  c^,  we  compare  the  oorres- 


vs.  CV/3K  is  available  in  Ref.  22.  The  results  are  plotted  in  I*ig.  3« 
Curves  (a) ,  (b) ,  and  (c)  are  the  results  of  the  present  calculations. 


(a)  using  0°K  elastic  constants  (15)  and  Houston*  s  approximation  to 


G-(oo)  ; 


(b)  using  Houston*  s  approximation  to  G-(u)) ,  adjusted  for  the  temperature 
variation  of  the  elastic  constants  by  (13)  with 


(25)  e 0/cq,  =  5.17  x  1011/oT,  and 


O51  =  (5-17  -  0.00266671)  x  1011  dynes/cm2  5 

(c)  using  the  AC  approximation  to  G-(w) ,  adjusted  for  the  temperature 


variation  of  the  elastic  constants  by  (18)  with 
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(26)  c0/ct  ~  x  IO^/oiji,  and 

cT  =.  (4.906  -  0.001787  T)  x  1011  dynes/em2. 

Curve  (f)  represents  a  smoothening  out  of  the  experimental  results. 
Points  marked  *x*  are  due  to  Bueken,  Clusius,  and  Woitinek2°  and  those 

OA 

marked  *o*  are  due  to  Meads,  Forsythe,  and  Giauque^  .  In  addition,  the 
calculated  results  obtained  by  Bhatia  and  Horton J  are  included  as 
curves  (d)  and  (e).  Both  were  obtained  by  using  Houston* s  approximation 
to  G(«o)  and 

(d)  0°K  elastic  constants  (20)  and, 

(e)  room  temperature  constants  (19). 

Theory  and  Experiment  Below  15°K 

Although  available,  no  experimental  results  have  been  shown  for 
the  temperature  region  below  15  °K.  The  work  of  Keesom  and  Kok25 
gives  a  sharp  maximum  in  the  G,.  -  T  curve  around  5°  followed  by  a 

minimum  around  5°,  These  findings  are  corroborated,  within  experimental 

26 

error,  by  the  results  of  Keesom  and  Pearlman  .  However,  these  two 

21 

groups  of  results  are  at  variance  with  those  of  Corak  et  al.  which 

g  o 

indicate  that  the  atomic  heat  obeys  the  T^  law  up  to  at  least  4  K ; 

consequently,  0^  is  constant  at  225.3°  to  4  °3£.  In  the  case  ot  Corak, 

it  must  be  noted  that  only  after  an  adjustment  of  the  1948  helium 

3 

vapour  pressure  -  temperature  scale  did  the  T  law  adequately  represent 
the  experimental  results.  However,  essentially  the  same  adjustments 
were  indicated  by  the  work,  along  quite  different  lines,  of  other 
investigators.  Thus,  Corak* s  corrections  ?;ould  appear  to  be  justified. 

A  small  part  of  the  discrepancy  between  theory  and  Corak* s  results 
may  have  its  origin  in  the  numerical  integration.  The  value  of 


-Sl¬ 


ob  tain  ed  by  numerical  integration  of  (1)  is  quite  sensitive  to  the 
frequency  interval  taken.  This  was  found  to  be  especially  so  for 
2  °K,  where  a  wide  discrepancy  between  9^  and  was  traced  to  this 
source*,  by  taking  smaller  integration  intervals,  an  improvement 
substantial  enough  to  give  agreement  between  these  two  values  was 
obtained. 

Theory  and  .ksgeriment  ^bove  15° 

A  feature  common  to  all  the  theoretical  curves  is  the  presence 
of  a  deep  minimum  around  20  °K.  This  minimum  is  primarily  due  to  the 
first  peak  in  the  G-(*>)  curve,  although  there  is  some  contribution  from 
the  peak  occurring  around  co(a^p)^'^  s  20  x  10^.  Thus,  because  both  of 
these  peaks  have  been  enhanced  by  the  AG  Method,  the  minimum  in  curve 
(c)  has  been  considerably  deepened.  On  the  other  hand,  it  is  largely 
because  of  the  second  peak  that  we  obtain  such  good  agreement  near 
60  °T,  for  curve  (c).  The  results  indicate  that  what  is  required  is 
a  redistribution  in  the  G-(w)  curve,  with  the  first  two  or  three  peaks 
somewhat  reduced  and,  the  contribution  of  the  portion  beyond 
-  20  x  IQ5  increased. 

Since  the  peaks  have  their  origin  in  the  approximations  through 
the  choice  of  directions  used,  it  is  quite  likely  that  the  actual 

is  flatter  than  thktj  has  been  used  here.  This  is  confirmed  by  the 
G(oj)  for  aluminum,  obtained  using  thousands  of  directions  .  There  are 
two  maxima-,  a  fairly  smooth  one  corresponding  to  the  transverse 
vibrations,  and  a  sharper  one  from  the  longitudinal  vibrations. 

Sane  of  the  discrepancy  between  theory  and  experiment  probably  has 
its  origin  in  the  elastic  constants  and  their  variation  with  temperature. 
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The  limitations  of  the  procedure  used  to  adjust  G(oo)  for  this  variation 
are  evident  from  the  assumptions  made  about  the  behaviour  of  the 
constants  with  temperature,  in  error  of  a  couple  of  percent  from 
this  source  alone  is  quite  conceivable. 

Another  factor  is  the  possible  influence  of  anharmonic  terras, 
whose  effect  would  naturally  be.  most  evident  at  higher  temperatures. 

We  note  that  curve  (b)  (for  which  G-(cc)  was  temperature  adjusted)  does 
not  remove  all  the  discrepancy  between  theory  (curve  a)  and  experiment 
at  the  high  temperature  end,  so  that  one  might  suspect  the  contribution 
of  anharmonic  tenas  to  the  specific  heat.  However,  since  the  temperature 
correction  itself  is  somewhat  uncertain,  we  should  not  jump  to  any 
conclusions  regarding  the  effect  of  anhaimonicity. 

Conclusions 

Taken  ov*er  the  whole  temperature  range,  the  Debye  characteristic 
temperatures  calculated  on  the  basis  of  the  Bhatia  theory  (curve  d)  and 
those  on  the  basis  of  the  Born  and  Begbie  theory  (curve  a),  using  the 
Houston  approximation  to  G(cj)  in  both  of  them,  appear  to  be  about 
equivalent.  The  former  theory  leads  to  slightly  better  results  at  the 
minimum  around  BO0,  but  beyond  about  80°,  the  latter  theory  fits  the 
experimental  results  more  closely,  but  the  differences  are  really 
quite  small. 

The  results  do  seem  to  indicate,  however,  that  a  three-term 
approximation  to  GW)  is  not  likely  to  lead  to  results  for  Q-p  that  are 
correct  r  the  whole  temperature  range.  The  number  of  directions  in 


such  approximations  is  too  small  to  give  adequately  the  frequency 
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distribution*  It  would  seem  likely  then  that  calculations  using  a 
larger  number  of  texms  in  G(co)  may  lead  to  a  G(c o)  raore  representative 
of  the  actual  picture,  and  thus  to  better  results  for  the  specific 
heats.  Such  a  calculation  might  indicate  significant  differences 
in  the  specific  heats  as  obtained  from  the  Bhatia  and  the  Born  and 
Begbie  theories. 


Accuracy  of  Work 

It  is  estimated  that  the  numerical  work  is  at  least  good  to  2$. 
Calculations  were  always  carried  out  to  five  figures,  while  all  graph 
readings  were  good  to  three  figures.  Intervals  for  any  numerical 
integration  were  always  sufficiently  close-spaced  to  capture  all  the 
details  of  the  variations  of  the  integrand. 
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APPENDIX  A*  Determination  of  Coefficients  of  Cubic  Equation  in  v^. 

Let  x,  y,  z  be  the  direction  cosines  of  the  wave-vector  k. 

=  xk;  k2  =  yk;  k^  =  zk 

In  the  long  wave  limit  (ka  0),  we  can  replace  cosines  and  sine3  of  ka 
by  their  Taylor  expansions  of  one  or  two  terms.  Thus  the  and  of  the 
secular  equation  (1.28)  become 

(1)  C12  =  1  -  (a 2k?lk)(xz  +  y2) 

C2^  =  1  -  (a2k2/4)  (y2  +  z2) 

C31  =  1  -  (a2k2/4)(z2  +  x2) 

=  xak//2;  S2  =  ya k/Jl;  S j  =  zak//2 

Substituting  these  expressions  into  the  D^  ,  we  obtain 

(2)  Dn  =  (k^/p^aCd  -  x2)  +  p(l  +  x2)} 

=  (k2/p)  +  x2(c1;l  - 

=  (k2/p)c^  {l  +  x2(ff+  i)  } 

B22  =  (k2/p)°44  I1  + 

=  (k Z/£)ci^  {l  +  z2(o*+  f ) 

D12  '  D21  = 

D23  =  D32  =  (k2/£^z^°44 

D31  =  D13  =  (k2/p)zx?c44 

wiisr©  0”  —  (o^  ^12  2o44^  / c44 

?  =  (oi2  +  tW|/c44 


-  «2 

D12 

D13 

>12 

D22  "  “2 

D23 

*13 

D23 

D33  -  ^ 

From 

(3) 


■ 


•I 


* 

'  •  1  :■ . :  ■  r  -  j  / 

x 

■ 

'  />-  l  <)!.,' . .  r*- 

(  '  :•  :  '  *'  -  ■  ■  ■  ..j 

:  -  : 


,•  '■  ' 
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we  obtain  the  cubic  equation  in  uj2 

(4)  J  +  Dzz  +  ^  +  a,2(DliD22  +  D22D33  +  ^  -  D§3  -  D^) 

^D11D22D33  +  2D12D23D31  ~  ■DX2I)33  ■  •D23D11  -  ■D31D22^  =  0 


Next,  lire  substitute  for  the  their  expressions  in  terras  of  the  elastic 
constants,  replace  or^/k^by  v^,  and 

X2y2  +  y2z2  +  Z2X2  g 

4x^y^z^  by  f 


We  obtain 

(5)  -  (Tif/p)(c11  +  2c^)  +  (vc i^/f>)2p  +  2(<r+£)  +  g  (cr+  2£)j 

-  (^/p)3!1  +  cr+  4  +  3^)  +  g  2£)^  =0 


which  we  write  as 

(6)  v6  ~  (v^/p)  (c^  +  20^)  +  (vci4^/p)2(p  +  tg)  -  +  sf  +  tg)  =  0 

where  p  =  3  +  2(cr+  £) 
r  =  <r  +  £  +  1 
s  =  o-2(<r+  3|)/4 
t  =  <r(cr+  2£) 

Therefore  M  =  (c^/p)2(p  +  tg) 

3ST  =  (c^/p)3(r  +  sf  +  tg) 

Finally, 

(7)  J(M/H)3/2cUl=  (p/c4^  Jl. 
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APPMDIX  B.  General  Depressions  for  (dnI)  /(dt^dffi) 

Writing  I  =  (P/q)°/2,  where 
P  -  a  -  g,  and 
Q  =  c  +  df  -  g, 

we  ha^e  for  the  general  expressions  for  (dnI)/(df%g^) : 


a  i 

3  P3/2d 

df  ’ 

2  <£>^ 

d  I  „ 

3P3/*  3P1/2 

a  g 

2  q5/2  2  Q3/2 

a2i 

15  p3//2a2 

df2 

'  4^7F- 

d2I 

15  pS/2d  9 

pi/% 

dfag 

'  4  q7/2  '  4 

q5/2 

<3  2I 

15  P3/2  9  P1/2  3  P-1/2 

dg2  ‘ 

4?72  zJl. 

2  4  Q3/ 2 

d3i 

105  ?3/2d3 

o,, 

%\ 

ii 

8  q972 

d  3I 

105  P3/2a2  ■ 

48  P1/,2d2 

a  f8dg 

8  q9/2 

8  Q7/2 

a3i 

105  P3//2d  , 

45  F17  2d 

9  p”^2d 

di'dg" 

8  Q9/2 

4^ 

8  q5/2 

d3i  . 

105  p3^2  135  : 

P1/2  .  27  P' 

■1/2  3  p-3/2 

ag3 " 

8  q?/2  8  i 

q7/2  '  8  q5/2  8  q3/2 

aii  . 

945  P3,/2d4 

at4  “ 

16 

a4i 

945  P3^2d3  ^ 

315  P^V 

5 

df°dg 

16  Q11/2 

16  q9/2 

$4I 

945  p3/2d3 

315  p-^a2 

4.5  P~1/2d2 

df~ag* 

2  16  q11/2 

8  q9  /2 

16  q7/2 

-57- 


a4i 

945  P3//2d 

945  P1//2d 

135 

p-i/2^  9 

p-3/% 

dido3 

16  q11/2 

16  Q9/2 

16 

q7/2  16 

Q5/2 

d4I  _ 

945  P3/2  315 

P1/2  45 

p-1/2 

9  P'3/2 

-S/H 

,  9  P 

dg4  ~ 

16  ft11'2  4 

q9/2  4 

^7T 

*  4  q3/2  * 

16  q3/2 

d5i 

10395  P3//2a5 

II 

I 

32  q13/2 

a5i 

10395  P3/2d4 

2855  P1/,2d4 

&f4dg 

32  q13/2 

32  q1;L/2 

d  5I  _  10395  +  3835  p1//8a5  315  p~-L-/2d0 

df3dgS  32  (-i13/2  16  :>11/  2  32  Qs/2 

(3  5I  _  10595  P5/8a3  _  8505  P^d2  ,  945  p~1//2d2  +  45  p~3/2d2 
dfSds3  32  o}'°//2  32  Ci11/2  32  Q,9/3  32 

d°I  .  _  10595  ?°//aa  t  2835  P^d  315  F~1,/3d  45  P~g//Sd  27  P-5^ 

dfdg4  32  Q,13/2  8  q11/2  32  q5/2 

<?5I  _  10395  ?3//3  14175  F1''3  t  315  p"1/2  45  P~S/2  (  135  p~3/2  |  45  P ~'?/2 

dg°  32  «L3/2  32  Q,11/2  4  cf/2  4  q7/2  32  Q5/2  32  Q3/2 
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APPENDIX  C 

Table  of 

Factors  Used  in 

Calculating 

I  1  c! Sl  bv  Ilethod 

(a)  Iron 

cos  0 

W 

X 

Y 

Z 

0.0 

1.33133 

-0.25000 

0.53479 

-0.25000 

0.1 

1.32143 

-0.24503 

0.52489 

-0.24020 

0.2 

1.29293 

-0.23040 

0.49639 

-0.21224 

0.3 

1.24943 

-0.20703 

0.45289 

-0.17033 

0.4 

1.19693 

-0.17640 

0.40039 

-0.12081 

0.5 

1.14383 

-0.14063 

0.34729 

-0.071381  j 

0.6 

1.10093 

-0.10240 

0.30439 

-0.029794 

0.7 

1.08143 

-0.065025 

0.28489 

-0.0022707 

0.8 

1,10093 

-0.03240 

0.30439 

0.0084396 

0.9 

1.17743 

-0.009085 

0.38089 

0. 0053728 

1.0 

1.33133 

0.000000 

0.53479 

0.0000000 

(b)  Silv  er 


cos  Q 

w 

X 

Y 

-  Z 

0.0 

1.03543 

-0.25000 

0.43797 

-0.25000 

0.1 

1.02553 

-0.24503 

0.42807 

-0.23979 

0.2 

0.99703 

-0.23040 

0.39957 

-0.21070 

0.3 

0.95353 

-0,20703 

0.35607 

-0.16719 

0.4 

0.90103 

-0.17640 

0.30357 

-0.11606 

0.5 

0.84793 

-0.14063 

0.25047 

-0.065469 

0.6 

0.80503 

-0.10240 

0.20757 

-0.023593 

0.7 

0.78553 

-0.065025 

0.18807 

0.0030900 

0.8 

0.80503 

-0.03240 

0.20757 

0.011928 

0.9 

0.88153 

-0.009025 

0.28407 

0.0066027 

1.0 

1.03543 

0.000000 

0.43797 

0.0000000 

i 
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( c )  Copper 


cos  8 

.w 

X 

Y 

Z 

0.0 

1.0306 

-0.25000 

0.42197 

-0.25000 

0.1 

1.0207 

-0.24503 

0.41207 

-0.23938 

0.2 

0.99219 

-0,23040 

0.38357 

-0.20916 

0.3 

0.94869 

-0.20703 

0.34007 

-0.16407 

0.4 

0.89619 

-0.17640 

0.28757 

-0.11134 

0.5 

0,84309 

-0.14063 

0.23447 

-0.059586 

0.6 

0.30019 

-0,10240 

0.19157 

-0.017424 

0.7 

0.78069 

-0.065025 

o 

o 

03 

D** 

r-t 

O 

0.008419 

0.8 

0.80019 

-0.03240 

0.19157 

,0.015398 

0.9 

0.87669 

-0.009025 

0.26807 

0.007826 

1.0 

:  1.0306 

0.000000 

0.42197 

0.000000 

(d)  Lead 


cos  9 

I 

X 

^ 7 

z 

0.0 

0.82397 

-0.25000 

0.35955 

-0.25000 

0.1 

0.81407 

-0.24503 

0.34965 

-0.23908 

0.2 

0.78557 

-0.23040 

0.32115 

-0.20803 

0.3 

0.74207 

-0.20703 

0.27765 

-0.16180 

0.4 

0.68957 

-0,17640 

0.22515 

-0.10790 

0.5 

0.63647 

-0.14063 

0.17205 

-0.055299 

0.6 

0.59357 

-0.10240 

0.12915 

-0.012931 

0.7 

0.57407 

-0.06503 

0.10965 

0.012303 

0.8 

0.59357 

-0.03240 

0.12915 

0.017925 

0.9 

0.67007 

-0.00903 

0.20565 

0.0087167 

1.0 

0.82397 

0. 00000 

0,35955 

0.0000000 
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APPENDIX  3Ka)Use  of  Ell in tic  Functions  to  Evaluate  J 


(a)  The  three  fundamental  elliptic  integrals  are: 

(l)  Normal  elliptic  integral  of  the  first  kind, 


dt 


■Gr. 


av 


O y  (1  -  t2)  (l  -  k^t2)  Jo  J~l  -  k2sin2V 

rai 

=  I  du  =  u^  =  sn  ~1(y,k)  =  3P0P,k) 
Jo 

(y  =  sin  f  ;  =  am  u^) 

(2)  Normal  elliptic  integral  of  the  second  kind, 


du 


Jjp-'l"  dt  =  jjl  -  k*sinI»'  d}/=  J  dneu 

5  E(n^)  E  S(am  u^,k)  S  E(*f>k) 

(3)  Normal  elliptic  integral  of  the  third  kind. 


f 


dt 


(l  -  of  2t2)/  (i  -»  t^)  (l  -  k2't'2T  J  (i  -  jl  -  k2sin2  r 


■I 


dV 


=  f  — — 5  J|(aBi  u1#0C2,k)  =TT(<PA2,k) 

1  4a  oC  sn  u 

(-s©  <  oc  2< 

When  y  =1,  the  integrals  are  said  to  he  complete.  In  that  case,  one  writes: 

-fT/g. 


^  f  =  f  du  =  F(tt/ 2,k)  s  K(k)  =  E 

4  7l  -  k^sin2V  J0 

(2)  j  x/l  -  k^sin2V  d*  =  |  dn^u  du  -  E(iT/2,k)  s  E(k)  ~  E 

(3)  f- _ ay ,  =  rK— & _ 

4  (1  -  c(2sin2y)Vl  -  k2sin2Y  1  1  -  o^srAi 


=  TT(-n/2,AZ,k)  =TJ<dZ,k) 
L  *'*•] 


.  i: 


oi 


:  O 


wwm 


J:  f  j> 
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APPSUDIX  D(b).  Uae  of  Elliptic  Functions  to  Evaluate  J 

(b)  Integration 


I.  n~  <  0 


•l  fe 


Wr) 

s.rr  xrl 


HZ 


dr 


■j. 


n/z 

.  _ 

L  ( i  -♦>  I  n*  J  stn'r)  3,z 


fTT  (2. 

I  nx  s.n*  V  dr 

Z  '  1  (I  <-  in‘|  ■m'-rl 


(See  p,  33-34  for  validity  of  this  expansion  of  numerator) 


=  f£  *  3n?-dL  ra 

Jln^l  Jo  (I  + 


c/  r' 

Inf  I 


■*'E  W 


-i'£  *ir$,  ■*'(*-  -S'V)- f  "‘-£'k 


--  -i'e  ♦ |  =£?  (e.  -  *"W  -  -l  <-£‘k 


w  here 

-T,  <  -  -A* 

II.  1  > 

n,1  ?  n*  >  o 

/TT/2 

T  - 

}  —  n,1  Sc/rtV 

J  "  K 

1  -  Hi 

(jr/Z 

df  j  1 

~  l 

/  -  k*  s^xrj  i 

I  1  n£| 


-TT/e 


|  -  n,  surt  r 

(  -  5unlr 


f  JTI Z  r  i  X  . 

=  «£  -  Zn,1  1  /  '  -  n-  dr 

7  dnl  JQ  J  i  -  n^surt  V 
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*  ?E  -  z"‘  £;l 
*  1  A:  J 


7TJ  2. 


d  V 


yp"-  n*  sun^y  )( *  -  1  r; 


7T  ^  2. 


sunly  ctr 


Jc  /(/  -  s^y)(  n^sw'r) 

=  <*£  -  Z n,2-  A  x  (aK)  +  Zn,1  J  x  Q  f  *  sn* - 

7  dnt  j  •  J0  t-  +  n*  $yxl 

=  a£  _  n,MK-e)_g»  ^  2n«  j9  9J  f 

*  -4*“  '  dry?  *  K  l  -  dEz  Jhzu 


where.  at  -  -  — Qi_ 
f-  nj 


-  ap  _  n,^11  ,  e)  2n1  d  n3  f  K_  _  TT  AoK-£>„ 

?  ^(K  }  2  ■  d*  »  U1-*1  t  TTuT^m^-*) 

-  aE  -  aif1  (K-  E)  +  2«/  A  [  Ai  -  2  A jiiii  ] 

7  '  •  an‘  L  h,1  e.  *■>,  n».  J 


?E 


^3(K-£>  +  ^3(k-£)  +  E  £>.<^-4) 


l  n» 


n-  » 


=£[**-  **«.<».*)] 


w 


h«re  4"=  ;  f- 


I-  Hxl 


V  =  s^'1 

n, 


TIL  I  ^  n*  7  n%x  ?  o 


r  ^  ,  r 

i  -  5^Ly 

r,z  wr* 

f  /  -  n,*"  5  m7*- 

)o  i  -  .«V  J 

i  ~  n*  swlr 

’  Jo  «-r 

/  -  nf  s^h  V 
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/  K 

-  rui^u  Au.  _  2.  n1  ~ 

*  Jo  •  dn?J 


ir/2. 


In 


«  a 


1  /f'  '  r\l  S  im  VTT*  '  ^  S  i^f'y) 

mz. 


dn 


y  (  r  -  n,1  3t^r)  0^  Jt/kxV 

=  %3  -  2m,2-  ^  ,  (a K)  +  zn*  £  a  f* _ »".*«■  ^ 

»  '  c><  *  Jc  I  -  n,1  1.  n,'  Sm'* 

=  _  "*  ^(e--&'Z'<)  2  «  J?  _J_  fK  s^V  dec 

4'1  ^  0<  ,-^y0  ,-^^u. 

w  h  *  r  e  oCr  =; 


n. 


=  4^-  -  x,tai(e--fe'\)  ?  ,  Ji  f  K  _  jr  />„(-»,  4) 

ttf*  f  ^  S^Jo^vI^T' 

•  iF*  -  z"'9e<U;j  -"Vo^d 

=  4^  -  n,V(  6-  4'  K)  ^  n^/e-^'K)  _  n^K  ^  tt  n,J  rt,,  (^-&> 
-ft  J^'z  nj 


If 

-  ".y( 

'-  $)l 

'  H-  4'V) 

>  -*'■*•*■  / 

-_  4\ 

-  l!ii  ( 

p-^ 

)  -  **  1 

4'1 

*  "  l 

)  Hj 

=  4J; 

(i“ 

n,lK  ) 
"*z  / 

-  tt  n,3 
t  r\3 

a  n? 


where  I-  -ft'^  rV*-..r^ 


I  -  rv 


a'.  _L_ 

*  I  -  n, 


i  ' 


=  S  isvij 


-  I 


Ao^-4)  -  %  [EF(V,4)  4.  KE(t-fe')  -  KP(y,4';] 


APPENDIX  E.  Id Jl  by  Method  III 
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1 


(b)  Copper 
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I  I  I 


(c)  Lead 
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(d)  Potassium 


■67- 


9- 
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CA 
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VO 

CA 

co 
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rH 
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CM 
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UN 

00 
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00 
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